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PROBLEM  DESCRIPTION 


In  the  event  of  a  war,  earthquake,  or  other 
similar  disaster  anywhere  in  the  world.  United  States 
military  transport  aircraft  may  respond  by  delivering 
equipment,  supplies,  or  personnel  to  the  affected  area. 
These  items  are  loaded  into  the  aircraft  taking  the 
physical  restrictions  of  maximum  weight  limit,  type 
and  size  of  equipment,  and  others. 

The  C-5  aircraft  is  the  largest  cargo  transport 
in  the  free  world  with  a  maximum  gross  weight  of 
760,000  pounds.  It  has  been  calculated  that  17  C-5s 
could  have  accomplished  the  entire  Berlin  Airlift. 

This  analysis  was  based  on  C-5A  characteristics  since 
it  is  likely  to  be  involved  in  any  large  scale  airlift 
operation. 

Another  physical  problem  comes  up  when  the 
distance  between  the  origin  base  and  the  destination 
base  is  too  far  for  the  aircraft  to  travel  without 
refuelling  at  a  location  in  between.  Bases  may  be 
located  along  the  way  or  it  may  be  more  efficient  to 
mid-air  refuel  using  special  aircraft.  The  tanker 
aircraft  flies  to  a  point,  meets  the  transport  aircraft, 
and  fuel  is  transferred  from  the  tanker  to  the  transport 
while  they  are  still  flying. 


The  objective  is  to  complete  the  move  of  the  cargo 
from  the  origin  to  the  destination  with  minimum  total 
fuel  consumption. 

Since  we  specified  the  use  of  the  C-5A  aircraft, 
we  need  to  determine: 

( 1 )  weight  of  the  cargo  load 

(2)  how  much  initial  fuel  to  load 

(3)  transport  route 

(4)  tanker  route 

(3)  and  (4)  are  determined  by  the  refuelling 

point. 

We  are  only  considering  the  single  refuelling 
case  where  there  are  two  legs  of  the  trip  for  the  cargo 
aircraft.  An  option  was  added  to  the  program  to  detect 
if  refuelling  was  not  necessary,  neither  mid-air  nor 


at  a  base 


ASSUMPTIONS 


(1)  The  total  cargo  can  be  divided  any  way  desired, 
i.e.,  we  can  ignore  the  fact  that  the  cargo  is 
loaded  on  pallets  and  that  we  can’t  split  pallets. 
Therefore,  the  cargo  weight  can  be  considered  a 
continuous  variable. 

(2)  No  size  restriction  about  the  load,  i.e.,  no 
part  of  the  cargo  is  too  big  to  fit.  This  is 
fairly  logical  for  the  C-5A  aircraft. 

(3)  The  only  limitation  on  the  weight  of  the  aircraft 
is  the  gross  weight  of  the  aircraft  at  take-off 
should  not  exceed  the  maximum  take-off  weight. 
Aircraft  Empty  Weight  +  Cargo  Weight  +  Fuel  Weight 

€  MTOW 

(4)  The  earth  is  a  perfect  sphere  with  radius  R  where 

R=39?0  status  miles 
=3404  nautical  miles 

(5)  Measure  of  "distance".  The  shortest  distance 
between  any  two  points  on  the  surface  of  a  sphere 
is  the  shorter  great  circle  arc  connecting  them. 
The  great  circle  is  the  circle  in  E^  which  has  a 
center  at  the  center  of  the  sphere  and  is  drawn 
on  the  surface  so  that  its  radius  is  equal  to  the 
radius  of  the  sphere. 


(6)  The  aircraft  follows  the  great  circle  arc  when 


4 


flying  from  one  point  to  another. 

(7)  Weather  is  assumed  to  be  negligible.  (Jet  streams 
or  storms  have  no  effect.) 

(8)  Aerial  refuelling  takes  a  negligible  amount  of 
time.  Therefore,  the  region  in  which  the  refuelling 
takes  place  is  considered  a  single  point. 


.v; 
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NOTATION/DEFINITIONS 

(1)  Only  C-5A  aircraft  are  available 

(2)  EW  s  transport  Empty  Weight 

(3)  MTOW  =  transport  Maximum  Take-Off  Weight 

(4)  =  transport  Maximum  Fuel  capacity 

(5)  H  «  tanker  Maximum  Fuel  capacity 

(6)  w  =  weight  of  the  cargo  to  be  transported 

(7)  gQ  =  initial  fuel  of  the  cargo  aircraft 

(8)  h  =  initial  fuel  of  the  tanker  aircraft 

(9)  R(gQ»w)  =  Range  of  the  transport  when  its  initial 

fuel  is  gQ  and  its  cargo  weight  is  w 

(10)  FC(go,wfd)  =  Fuel  Consumed  by  the  transport  when 

it  flies  a  distance  d  and  its  initial 
fuel  is  gQ  and  its  cargo  weight  is  w 
(note:  d  must  be  «£  R(gQ*w)  ) 

(11)  FN(w,d)  =  exact  amount  of  Fuel  Needed  to  fly  a 

distance  d  when  the  cargo  weighs  w 
(note:  d  must  be  £  ^ 

(12)  L(hQ)  =  range  of  the  tanker  when  its  initial  fuel 

is  ho 

(13)  FCT(h0,d)  =  Fuel  Consumed  by  the  Tanker  when  it 

flies  a  distance  d  and  its  initial 
fuel  is  li 

o 

(note:  d  must  be  ^  L(^0)  ) 


(14)  FNT(d)  =  exact  amount  of  Fuel  Needed  by  the  Tanker 

to  fly  a  distance  d 

(note:  d  must  be  ^  L(H  v)  ) 

(15)  DQd  =  the  shortest  distance  between  the  origin  base 

and  the  destination  base  measured  along  the 
arc  of  the  great  circle  connecting  the  two 
points. 

(16)  (0Q»0O)  =  the  coordinates  of  the  origin  base  in 

spherical  coordinates  (longitude  and 
latitude) 

(ed,0d)  =  the  coordinates  of  the  destination 
(©b»0b)  =  the  coordinates  of  the  tanker  base 
(9,0)  =  any  point  on  a  sphere  (earth) 

(17)  d-j  =  d(0o,0o,©,0)  =  the  distance  between  the  origin 

and  the  point  (6,0) 

d2  =  (0d»0d,e’^ 

d5  =  (eb,0b,e,0) 

where  d(0.  ,0.  ,9  .,0  .)  =  the  shortest  distance 

function  between  the  two  points  (9^,0^)  and 

(9.,0,)  measured  along  the  arc  of  the  great 
J  0 

circle  connecting  the  two  points. 

(18)  GW  =  Gross  Weight  of  the  transport 

=  EW  +  +  w 

°o 

GW^ax  -  Maximum  Gross  Weight  of  the  transport 

(19)  g  .  minimum  (Pmax,OV(max-w) 


MEASURE  OP  DISTANCE 


The  closest  measure  of  distance  models  the  earth 
as  a  sphere.  While  it  is  a  spheroid,  and  not  a  perfect 
sphere,  it  is  close  enough.  One  travels  on  a  curve 
when  going  from  one  point  to  another  on  the  surface 
of  the  sphere.  The  curve  which  has  the  shortest 
distance  between  two  points  on  the  surface  of  the 
sphere  is  the  great  circle.  The  great  circle  is  the 
circle  in  E  which  has  a  center  at  the  center  of  the 
sphere  and  is  drawn  on  the  surface  so  that  its  radius 
is  equal  to  the  radius  of  the  sphere.  Therefore,  the 
shortest  distance  between  two  points  on  a  sphere  is 
the  length  of  the  shorter  great  circle  arc  connecting 
them. 

Several  methods  have  been  reported  in  the 
literature  to  calculate  this  distance.  All  of  the 
methods  use  two  facts: 

(1)  the  great  circle  distance  between  two  points 
on  a  sphere  is  directly  proportional  with 
the  radius  of  the  sphere. 

(2)  the  distance  between  two  points  on  a  unit 


sphere  is  identical  to  the  angle  (measured 
in  Radians)  between  the  two  normals  of  the 
sphere  at  these  points. 


So  the  distance  =  R** 
where  R  =  the  radius  of  the  sphere,  and 

=  the  angle  between  the  two  normals  at 
these  two  points  (measured  in  Radians), 
The  methods  differed  in  the  way  they  found  the 
angle  ©<. 

By  representing  X  and  Y  in  spherical  coordinates 
as  (R, 9-1,0.,)  and  (R,0?l0?),  respectively,  instead  of 
cartesian  coordinates  (x1tx2,x^)  and  (y -,,y 2'^^ »  the 


X1 

= 

R 

sin01 cos©, 

yi  = 

R 

sin0?cos9 

x? 

= 

R 

sin0, sin9, 

y?  - 

R 

sin0?sin9 

x3 

= 

R 

COS0-, 

ii 

R 

cos02 

11 

* 

11 

=  the  Euclidean 

norm 

11 

X 

11 

=  11  Y  11  =  R 

pAC 

X  .  Y 

) 

11  X  11  11  Y 

11 

11 

X 

11 

*  11  Y  11  = 

R2 

X  *  Y  =  R?( sin01 sin0pCos9, cos©?  +  sin0, sin02Sin9, sin9 
+  COS0,COS0p) 

=  R2(sin01 sin0p( cos9, cos02  +  sin9,sin9p) 


+  cos0.,cos09) 


If  we  recognize  that: 


cos(02  -  9-j )  *  cosS^osQp  +  sin©  .|  sin©  ? 

We  get: 

X  *  Y  =  R2(sin01sin02cos(©2  -  0^)  +  cos0.jcos02) 
cos  =  3in01sin02cos(@2  -  9^)  +  cos0.jcos02 

where  0  =  the  geographical  longitude 

0  =  PI/2  -  the  latitude  =  the  colatitude 

To  take  advantage  of  computer  routines,  we  used 
the  following  version  of  this  measure: 


J  1  -  (sin01sin09cos(01-9o)  +  cos01cos0o) 

d  =  R  *  arctan(  - ' - - - - — - - - - 

sin01sin02Cos(01-02)  +  0030^0302 


where  9  =  longitude  *  PI/180 

0  =  PI/2  -  latitude  *  PI/ 180 
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FUEL  CONSUMPTION  FUNCTIONS 

This  section  develops  the  relationships  between  cargo 
weight,  fuel  consumption,  distance  travelled,  and  initial 
fuel. 

Charts  are  available  from  the  U.S.  Air  Force  for 
aircraft  performance  at  different  gross  weights,  speeds, 
and  altitudes.  These  charts  need  to  be  transformed  into 
mathematical  formulae  before  they  can  be  used  in  a  model. 
An  example  of  this  data  is  shown  in  Figure  t.  Different 
curves  are  given  for  different  weights  and  show  the 
distance  travelled  per  1,000  lbs  of  fuel  used  at  a  given 
speed  and  at  the  particular  gross  weight.  It  is  assumed 
that  the  transport  will  operate  at  the  99%  maximum 
specific  range  to  maximize  the  distance  travelled  per 
unit  of  fuel  burned.  The  distance  in  miles  per  1,000  lbs 
of  fuel  burned  is  plotted  against  gross  weight  GW. 

The  y-coordinate  is  called  MPF(GW)  because  it  is  a 
function  of  the  gross  weight.  Two  kinds  of  curves  to 
fit  those  points  are  tried:  linear  and  quadratic. 

Both  fit  well  as  seen  in  Figure  2. 
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For  the  linear  fit  we  have: 

MPF(GW)  =  aQ  +  a.,GW 
where : 

ao  =  36.2829 
a1  =  -0.0270 

and  q  =  correlation  coefficient  =  -0.9918 
For  the  quadratic  fit  we  have: 

MPF(GW)  =  bQ  +  b.jGW  +  b2GW? 
where : 

bQ  =  43.7616 
b1  =  -0.0576 
b?  =  2.94  x  10"5 
and  q  =  -0.9983 

We  will  use  the  linear  fit  to  find  the  distance 
the  transport  oan  travel  with  initial  fuel  gQ  and 
cargo  weight  wQ.  We  can  also  use  it  to  determine  how 
much  fuel  is  used  to  travel  a  distance  d  with  cargo 
weight  wQ  and  initial  fuel  gQ. 

Note  that  at  any  time  in  flight: 

GW  =  EW  +  wQ  +  f 

where  f  =  the  instantaneous  amount  of  fuel, 
f  changes  due  to  burning  fuel  in  flight;  therefore, 
GW  changes  too.  EW  and  wQ  remain  constant. 
Therefore, 


dGW  =  df 
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MPF(GW)  *  miles  per  1000  lb  of  fuel 
burned  for  a  given  gross 
weight 
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GW  =  Gross  Weight  in  1000  lb 


Figure  1.  MPF  vs.  GW 


U 


Let  R(g0»*0)  =  the  range  of  the  transport  when  its 

initial  fuel  is  gQ  and  its  cargo  weight 

is  w  . 

EW+wrt +grt 

then  R(gn,wn)  °  0  MPF(GW)dGW 

0  0  EW+w. 

o 

MPF(GW)df 

MPF(EW  +  wQ  +  f )df 
Using  the  linear  fit  for  MPF, 

^1 

R(g0.w0)  =  (aQ  +  ai  wo  +  5  g0)go 
where  a^  *  aQ  +  a-^EW 

Now  we  must  find  the  fuel  consumed  when  the  aircraft 
flies  a  distance  d  (d*R(gQ,wo))  with  initial  fuel  go 
and  cargo  weight  wQ. 

Let  g^  =  final  amount  of  fuel  left,  and 

FC(g0,wQ,d)  =  fuel  consumed  when  initial  fuel  is 

gQ,  cargo  weight  is  wQ,  and  the  distance 
flown  is  d. 
then  FC(g0,w0,d)  =  gQ  -  gf 


g. 


0 


Therefore,  to  find  PC  we  need  to  find  g^,  which  can 
be  done  by  solving  for  gf  in  the  following  equation: 

g 

d  »  /  0  MPF(EW  +  wQ  +  f )df 
gf 

For  linear  fit  of  MPF  we  have: 

d  »  f0  (a  +  a1  (EW  +  f  +  w  ))df 
gf 

=  a0(g0  “  gf>  +  a1(EW  +  wo>*go  -  gf>  +  '1(«o  -  4 

0  =  “1gf  +  <ao  +  al(EW  +  %))gf 

-  <aos0  *  ai(EW  *  vo)go  *  -  d 


-  gf  +  aSf  -  d '  *  0  where 

This  equation  is  quadratic  in  gf. 
to  this  equation: 
a 


gf  =  -• 


+  V 


a  +  2a-jd’ 


l1 


a  =  aQ  +  a-j  (EW  +  wo) 

To  find  the  solutions 


(a  +  a^  V 


2a  .d 


a 


Only  the  positive  solution  makes  physically  possible 
so: 


Since  the  fuel  consumed  =  FC(go,wQfd)  =  gQ  -g^,  then 

FC(g0,w0,d)  =  g0  ♦  -  - 

a1 

where  a  =  a  +a1(EW+w) 
o  i  o 

Let  FN(wQ,d)  =  the  exact  amount  of  fuel  needed  to  fly 
a  distance  d  when  the  cargo  weight  is  wQ.  We  need  to 
find  FN.  For  a  linear  fit  of  MPF  we  have: 

a1 

R(g  ,w  )  =  d  =  (a'  +  a^w  +  -  g  )g 
VBo*  o'  v  o  1  o  2 

With  FN  =  gQ  and  solving  R  for  gQ  we  get: 


(a  +  a1gQ)?  -  2a^d 
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MODEL  AND  PROPOSED  SOLUTION  METHODS 

The  case  where  there  is  one  aircraft  with 
predetermined  cargo  weight  and  initial  fuel,  but 
needs  to  be  refuelled  is  a  subproblem  of  the  much 
larger  problem  involving  many  aircraft  and  multiple 
refuelling.  It  is  one  of  the  building  blocks  to  the 
solution  of  the  larger  problem. 

This  problem  can  be  broken  down  into  two 
subproblems  itself.  The  first  subproblem  can  be 
defined  as  finding  the  initial  fuel  needed  and  other 
fuel  needed  to  complete  the  mission  given  a  feasible 
refuelling  point.  The  second  subproblem  is  to  find 
the  location  of  a  refuelling  point  such  that  the  total 
fuel  used  is  minimized.  Using  the  results  from  the 
first  subproblem  as  input  for  the  second  subproblem, 
and  using  the  results  from  subproblem  two  as  input 
for  subproblem  one  in  turn,  an  optimal  solution  can 
be  found  which  does  not  require  us  to  predetermine 
the  initial  fuel  amount. 

Given:  EW,  MTOW,  GW,  and  FM  from  the  aircraft 

specifications  and  wo  from  the  mission 
requirements  we  can  determine  the  allowable 
fuel  loads  for  take-off  and  mid-air  refuelling. 


Let  FO  =  Fuel  weight  possible  to  add  to  aircraft  at 
take-off 

=  MTOW  -  EW  -  wQ 

GO  =  Maximum  fuel  weight  possible  at  mid-air 

refuelling 

=  GW  -  EW  -  w„ 
o 

GM  =  Maximum  fuel  weight  capacity  at  take-off 
=  min(FO,FM) 

GC  a  Maximum  fuel  weight  capacity  at  mid-air 
refuelling 
=  min(GO,FM) 

Given  a  suggested  refuelling  point  (9,0),  the  distance 
to  the  respective  bases  can  be  determined  with: 

LI  a  distance  to  the  origin  base 
D2  a  distance  to  the  destination  base 
D3  =  distance  to  the  tanker  base 
The  ranges  from  the  bases  may  be  determined  by: 

R(g0tw0)  =  R(GM,wq)  a  Range  from  origin  base 

=  R(GC,wq)  a  Range  from  destination  base 

a  R(Hmax>w0)  =  Range  from  tanker  base 

- - - 

(1/?  because  of  return 
trip) 

Now  the  feasibility  of  the  refuelling  point  can  be 
determined  by  comparing  the  distances  to  the  ranges. 
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For  feasibility: 

D1  ^  Rl,  D2  ^  R2,  D3  ^  R3#  and  the  distance  from 
the  origin  base  to  the  destination  base  =  Rl  +  R? 
Given  a  feasible  initial  point  we  may  find  the  minimum 
combination  of  initial  fuel,  transferred  fuel,  and  tanker 
fuel. 

Min(Total  Fuel  Cost  for  Subproblem  1 )  =  gQ  +  hQ  = 

J  (C  -  J  (a  +  algo)^  -  2a., d7  )*  +  2a "d3 


where  C  =  \/  a^  +  ?a1D2  +  V  a^  +  2a^D3 

We  search  for  the  minimum  total  cost  over  the  range 
from  Low  =  FN(wo,D1)  to  High  =  GM.  A  Golden 
Section  Search  is  used  in  the  program  to  find  this 
minimum.  The  accuracy  is  set  at  0.00010.  The  resulting 
Total  Fuel  Cost  is  checked  against  previous  solutions 
for  change.  If  there  is  little  change,  then  the 
problem  is  solved.  The  accuracy  is  specified  in  the 
initial  conditions.  A  large  TFC  has  been  defined  to 
caus  the  problem  to  cycle  at  least  once.  If  the  change 
is  significant,  then  go  to  subproblem  two. 


Subproblem  two  involves  finding  the  location  that 
minimizes  the  sum  of  fuel  costs  for  the  four  legs  of 
the  mission.  The  four  fuel  costs  are: 

FN(0,D3)  =  Fuel  needed  by  the  tanker  to  return 
to  base 

FN(wq,D2)  =  Fuel  needed  by  the  transport  to  get 

to  the  destination  from  the  refuelling 
point 

FN(Wp,D1)  =  Fuel  needed  by  the  transport  to  get 
to  refuelling  point  with  cargo  and 
unused  fuel. 

where  Wp  =  wQ  +  (gQ  -  FC(go> wq,D1 )) 

a' 

and  FC(go,wo,D1)  =  gQ  +  wQ  +  - 

a1 

FN(w,D3)  =  Fuel  needed  by  the  tanker  to  get  to 

the  refuelling  point  with  net  fuel  needed 
to  transfer 

where 

w  =  FN(0,D3)  +  FN(wo,D?)  -  (gQ  -  FC( gQ, wQ, D1 ) ) 

Subproblem  two  is  subject  to  the  same  feasibility 
constraints  as  the  original  problem. 


Finding  the  optimum  location  involves  searching 
along  0  and  0  in  an  iterative  approach  to  find  the 
minimum  fuel  cost.  Limits  for  the  different  functions 
were  hard  to  find  so  the  initial  feasible  point  is 
used  as  a  beginning.  0  and  0  are  changed  in  a 
systematic  manner  and  whichever  yields  the  lower  fuel 
cost  is  selected. 

(0,0)  yields  a  cost  F2 

(0+0, 0+d)  yields  a  cost  F2* 

If  (0+0, 0+0)  not  feasible,  try  another  point 

If  F2*  >  F2,  try  another  point 

If  F2 ’  <  F2,  F2  =  F2’  and  0  =  0+0,  0  =  0+0 

and  try  another  point  in  the  same 
direction 

If  F2*  =  F2,  then  reduce  0  and  try  another 
point  in  same  direction 

Similarly,  when  conditions  lead  us  to  try  another 
point: 

(0-0, 0+0),  (0-0, 0-0),  and  (0+0, 0-0) 

In  each  case,  0  returns  to  the  original  value  for 
the  new  direction.  After  one  complete  pass,  0  is 
checked  to  see  if  it  is  below  the  specified  accuracy 
limit.  If  it  is,  then  the  second  subproblem  is  solved. 
If  not,  then  reduce  the  original  value  of  0  and  make 
another  pass.  In  each  direction,  when  ties  occur, 

0  will  reduce  to  the  specified  accuracy  before  trying 


PROGRAM 


The  main  body  of  the  program  consists  of  defining 
of  constants;  input  of  locations,  cargo  weight,  and 
accuracy  requirements;  gosubs  for  the  subproblems; 
and  various  output  messages. 

The  transport  aircraft  and  tanker  aircraft  are 
assumed  to  have  identical  performance  charts, 
fuel  Cost  and  Fuel  Needed  subroutines  for  the  tanker 
would  have  to  be  changed  if  this  were  not  true. 

The  constant  C  would  also  have  to  be  changed. 

The  longitudinal  difference  between  any  two 
points  is  restricted  to  180  degrees  and  the  latitude 
is  restricted  to  North  latitudes.  Longitude  and  latitud 
are  converted  to  radians  using  the  formulas: 

9  =  (longitude  diff erence)*Pl/180  and 
0  =  PI/2  -  latitude*PI/180.  The  longitude  of  the 
origin  base  is  always  transformed  to  0. 

A  map  of  the  U.S.  and  Europe  is  available  at  the 
end  of  the  program  which  display  coordinates  from 
108°W  to  48° E  and  0°N  to  70°N.  This  will  allow  display 
of  most  U.S.  to  Europe  or  Middle  East  missions  and 
their  return  flights. 

The  longitude  of  the  tanker  base  should  be  in  the 
same  direction  from  the  origin  base  as  the  longitude 
of  the  destination  base. 
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The  following  lines  may  be  hard  to  read  in  the 
program  since  the  printer  also  printed  on  the  perfor¬ 
ations  separating  the  pages: 


1895  REM  CHECK  FOR  CORRECT  INPUT 
3200  GOSUB  04400  :  REM  DISTANCE(TA, PA) 

6400  GR=GM  :  WR=WO 

10400  G0=XB  :  SF=(A+A1*XB)f2-2*A1*D1 

15300  IF  F2=FF  THEN  TY=1 

22620  REM  WHICH  IS  SLIGHTLY  DIFFERENT  FROM  SUBPROBLEM 
ONE  BECAUSE  OF 

28030  REM  THIS  ROUTINE  MUST  BE  EXPANDED  TO  ALLOW 
ROUTES  ACROSS  THE 

35200  IF  (SJ$= ”E")  AND  (SL$=HW")  AND  (SN$="E")  AND 
(SQ$= "E")  THEN  L6=L1-CA 

60533  POKE 1 99 1 , 4  3 : POKE2000 , 4 3 : POKE20 1 2 , 4  3 : POKE202 1 , 4  3 
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100  EU-370.0  QQ-0  :  Al— 0. 026785  AO-36.1459  :  MT-790.0 

200  014=760.0  HM-300.0  FM-300.0 

300  AP=A0+A1*EW  PI-3.14139  :  R-3960. *5280. /6080.  ;  K-0  F0-1000000. 

310  REM  EW-EMPTV  WE  I  OUT/ 1000  LBS  MT-MAXIMUM  TAKEOFF  ME1OHT/1000  LBS 

320  REM  '3W-MAXIMUM  (3R08S  WEIGHT/ 1000  LBS  FM-MAX1MUM  FUEL  CAPRCITV/1000  LB": 

3 30  REM  HM-MAX I MUM  FUEL  CAPACITY/ 1000  LBS  OF  TANKER 

335  REM  FM-MAXIMUM  FUEL  CAPACITY/ 1000  LBS  OF  TRANSPORT 

340  REM  A 1 -LI NEAR  SLOPE  FROM  ENGINE  PERFORMANCE  CURVES 

350  REM  AO-LINEARTR  INTERCEPT  FROM  ENGINE  PERFORMANCE  CURVES 

368  REM  QQ=INITIAL  CASE  SETTING  FOR  POSSIBLE  LOCATION  SCHEMES 

370  REM  AP-FUEL  REQUIRED  TO  MOVE  EMPTY  AIRPLANE  R-8PHERICAL  CONSTANT 

380  REM  K- ITERATION  COUNTER  F0-INITIAL  TOTAL  COST  SET  VERY  LARGE 

408  PRINT  "3" 

410  PRINT  “THIS  PROGRAM  WILL  NOT  WORK  FOR" 

415  PRINT  "SOUTH  LATITUDES  OR  FLIGHT  PLANS" 

420  PRINT  "CROSSING  THE  INTERNATIONAL  DATE  LINE." 

425  PRINT  "FUTURE  MODIFICATIONS  WILL  CORRECT" 

430  PRINT  "THIS  SITUATION."  '  PRINT  PRINT 
440  PRINT  "ENTER  COORDINATES" 

450  PRINT  "IN  DEGREES ,  THEN  DIRECTION" 

455  PRINT  "EXAMPLE:  20  C RETURN)  E  (RETURN) " 

460  PRINT  "FOR  A  LONGITUDE  OF  20  DEOREES  EAST” 

470  PRINT  PRINT  :  PRINT 

500  INPUT  "ORIGIN  LONGITUDE", LI  INPUT  "DIRECTION  E/W" ; SJB  PRINT 
550  INPUT  "ORIGIN  LATITUDE" , HI  :  INPUT  "DIRECTION  N/S" ;8K* : PRINT 
600  INPUT  "DESTINATION  LONGITUDE". L2  INPUT  "DIRECTION  E/W" ;SL*: PRINT 
650  INPUT  "DESTINATION  LAT I  TUBE ” ; M2  INPUT  "DIRECTION  N/S" i SM# : PRINT 

700  INPUT  "TANKER  BASE  LONGITUDE", L3  =  INPUT  “DIRECTION  E/W" ;SN»: PRINT 
750  INPUT  "TANKER  BASE  LAT I  TUBE " ; M3  INPUT  "DIRECTION  N/S", SP*  PRINT 
800  INPUT  "CARGO  WEIGHT  (IN  POUNDS>“,WO 
900  PRINT  "FUEL  COST  ACCURACY! TO  NEAREST  POUND  =  1>": 


CTO  NEAREST  10  POUNDS" 
10V';E1 


910  PRINT 
920  INPUT 
930  PRINT 

1000  PRINT  "LOCATION  ACCURACY  DESIRED" 

1010  PRINT  "(TO  NEAREST  DECREE  -  I  >" 

1015  PRINT  "CTO  NEAREST  MINUTE  -  2  V 

1020  PRINT  "CTO  NEAREST  SECOND  -  3  >" 

1030  INPUT  "SELECT  ACCURACY", RZ  PRINT 

1040  WO-WO/1000 

1045  E1=E1/1000 

1050  IF  RZ-1  THEN  E2-1PPI/180 

1055  IF  RZ-2  THEN  E2-1/60BP 1/180 

1060  IF  RZ»3  THEN  E2»1/360*P 1/100 

1070  PEM  CONVERT  INPUTS  TO  LIKE  UNITS  FOP  CALCULATIONS 

1100  FO-MT-EW-WO 

1200  O0-OW-EW-WO 

1300  IF  FMOFQ  THEN  GM-FM 

1350  IF  FM>FO  THEN  OM-FO 

1400  IF  FMC-00  THEN  GC-FM 

1450  IF  FM;60  THEN  GC=O0 

1500  H-AF+R1PWO 

1600  PEM  SET  LIMITS  ON  FUEL  TO  DETERMINE  FEASIBILITY 
1700  INPUT  "PROPOSED  LONG I TUDE " ; CA  :  INPUT  "DIRECTION  E/W" , SO* : PRINT 
1800  INPUT  "PROPOSED  LATITUDE";CB  INPUT  “DIRECTION  N/S" , SRB  PRINT 
PRINT 

REM  CONVERT  TO  RADIANS 


1810  PRINT 
1850  GOSUB  28000 
1855  py«0 

I860  IF  ( QQ< 1 >  OR  <QQ>12)  THEN  PV-1 
1870  IF  8KB- "S 11  OR  SM*-"S"  OR  SP*-"S"  OR  SR*«"S"  THEN  PV-1 
1880  IF  PV-1  THEN  PRINT  "CHECK  INPUTS 
1885  IF  PV-1  THEN  PRINT  :  PRINT 
1890  IF  PV-1  THEN  OOTO  410 
1A95  PFM  rucrif  cno  moBCfT  luaiiT 


PARAMETERS  VIOLATED" 


1900  GOSUB  03000  REM  CHECK  FOR  FEFlSIBILITV 

1950  IF  IN-45  THEN  GOTO  410 

2000  IF  INO40  AND  INO50  THEN  GOTO  01700 

2005  REM  IF  INITIAL  LOCATION  IS  NOT  FEASIBLE-TRV  AGAIN 

2010  IF  IN-40  THEN  GOSUB  18000 

2015  REM  SOLVE  SPECIAL  CASE  UHEN  REFUELLING  NOT  HEEDED 

2020  IF  IN-40  THEN  GOTO  2775 

2025  REM  PRINT  ANSWER  TO  SPECIAL  CASE 

2030  REM 

2035  REM 

2100  GOSUB  07500  REM  SOLVE  SUBPROBLEM  ONE  -FOR  INITIAL  FUEL  AND  TRANSFER  FUE. 
2110  REM 
2115  REM 

2208  Z=F-F0  Zl-ABS(Z)  IF  ZK-E1  THEN  GOTO  02700 

2210  REM  IF  NEW  SOLUTION  WITHIN  ACCURACV  LIMITS  THEN  PRINT  SOLUTION 

2300  F0-F 

2305  PRINT  :  PRINT 

2310  PRINT  "TRIAL  SOLUTION  TO  SUBPROBLEM  ONE" 

2320  PRINT  "INITIAL  FUEL" ; 00*1000 
2330  PRINT  "TANKER/TRANSFER" ; HO*1000 
2340  PRINT  " TOTAL "IF* 1000: 

2350  PRINT:PRINT 

2360  PRINT  "WORKING  ON  PROBLEM  TWO": 

2370  PRINT 

2330  REM  PRINT  SOMETHING  SO  OPERATOR  KNOWS  PROGRAM  IS  UORKINO 
2385  REM 
2390  REM 

2400  GOSUB  11700  REM  SOLVE  SUBPROBLEM  TWO  TO  FIND  BEST  LOCATION 
2410  REM 
2500  K-K+l 

2510  REM  UPDATE  ITERATION  COUNTER 

2520  REM 

2530  REM 

2600  GOTO  02100 

2610  REM  RETURN  TO  SOLVE  SUBPROBLEM  ONE  FOR  NEW  LOCATION 
2620  REM 
2630  REM 

2760  GOSUB  50500  :  REM  CONVERT  ANSWER  TO  LONGITUDE  AND  LATITUDE 
2750  PRINT  "3"  PRINT  "MINIMUM  COST  SOLUTION" 

2775  PRINT  "OPTIMUM  REFUELING  POINT”  PRINT  L9; D7*,H8, D0* 

2730  PRINT  "INITIAL  FUEL ", 00* 1 000 

2735  PRINT  "TRANSFER  FUEL", <F-G0-F7-F3)*1000 

2787  PRINT  "TANKER  FUEL" , <F7kF3>*1000 : 

2790  PRINT  "TOTAL  FUEL  COST",F*1000 
2300  PRINT  " ITERATIONS" JK 

2310  PRINT  "ORIOIN"  PRINT  LI ; SJ*,Hl ; 3K*  PRINT  "DESTINATION- 
2815  PRINT  L2;SL*,H2;SM* 

2320  PRINT  "TANKER  BASE"  PRINT  L3;SN*,H3,SP* 

2830  PRINT  "CAROO  WEIGHT", WO* 1000, "LBS" 

2840  PRINT  PRINT  PRINT 

2845  REM  PRINT  FINAL  ANSWER  AND  REPEAT  ORIGINAL  INFORMATION 

2830  PRINT  "A  MAP  OF  U.S.  TO  EUROPE  AND  WEST  ASIA- 

2881  PRINT  "IS  AVAILABLE.  DO  VOU  WISH  TO  SEE'’"  INPUT  ”^V^N)";MP* 

2835  IF  MP*-"V"  THEN  GOSUB  60500 

2890  REM  MAP  DEMONSTRATION  AVAILABLE 

2895  GOTO  2750 

2900  STOP 

2910  REM 

2920  REM 

2930  REM 

2940  REM 

3000  REM  FEASIBLE 

3010  REM  THIS  ROUTINE  DETERMINES  IF  THE  PROBLEM  IS  FEASIBLE  FOR  THE  GIVEN 
3020  REM  PARAMETERS, 

3100  TA-TH  ■  PA-PH 

~  nopiip  .  pcm  tuSMNCElXA  pa  ,  _ _ _  — 


3210  REM  CALCULATE  THE  DISTANCE  FROM  THE  TRIAL  POINT  TO  THE  THREE  BASE  PC 
3300  GOSUB  06300  REM  RANGE 

3310  REM  CALCULATE  THE  RANGE  OF  AIRCRAFT  FROM  THE  THREE  BASE  POINTS 
3320  WL-WO  DL-01  OOSUB  08200  :  REM  FNTR(UO.Dl) 

3375  IN-0 

3400  IF  c  D 1 >R 1 >  OR  <FT>QM)  THEN  IN-10 

3500  IF  D2>R2  THEN  IN-20 

3680  IF  D3>R3  THEN  IN-30 

3610  IF  09<R1  THEN  IN-40 

3620  IF  09>R1+R2  THEN  IN-45 

3700  IF  IN-6  THEN  IN-50 

3800  IF  IN-10  THEN  r*»"TOO  FAR  FROM  ORIGIN"  ' 

3900  IF  1N=20  THEN  I*=”TOO  FAR  FROM  DESTINATION” 

4000  IF  IN-30  THEN  IS- "TOO  FAR  FROM  TANKER  BASE” 

4010  IF  IN-40  THEN  If»”REFUEL INO  NOT  NEEDED” 

4020  IF  IN-45  THEN  IS-” INFEASIBLE-ORIGIN  TO  DESTINATION  TOO  FAR" 

4100  IF  IN-50  THEN  IS--INITIAL  LOCATION  IS  FEASIBLE” 

4200  PRINT  IS  PRINT  PRINT 

4300  RETURN 

4310  REM 

4320  REM 

4330  REM 

4400  REM  DISTANCE<TA,PA> 

4410  REM  DISTANCE  CALCULATIONS 
4500  TC-T0  :  PC-PO 

4510  IF  ABS<TC-TA)>PI/2  THEN  OOSUB  19000  REM  SPECIAL  CASE 

4520  IF  RBS<TC-TA»PI/2  THEN  GOTO  4700 

4600  OOSUB  03300  =  REM  CALC-LENQTW  ALCTHI , PHI , TH2. PH2> 

4700  Dl-AL  TC-TD  ■  PC-PD 

4710  IF  ABS<TC-TA»PI72  THEN  OOSUB  19000  REM  SPECIAL  CASE 
4720  IF  ABS<TC-TA»PI72  THEN  OOTO  4900 
4800  GOSUB  03300  :  REM  CALC-LENOTH 
4900  D2-AL  ■  TC-TB  :  PC-PB 

4910  IF  ABS(TC-TAJ>PI/2  THEN  OOSUB  19000  REM  SPECIAL  CASE 
4920  IF  ABS<TC-TA»PI/2  THEN  OOTO  3300 
5200  OOSUB  05300  :  REM  CALC-LENOTH 
5300  D3-AL 

5310  TA-T0  PA-PO  TC-TD  PC-PD 

5315  IF  ABS(TC-TA)>PI/2  THEN  OOSUB  19000  REM  SPECIAL  CASE 

5320  IF  ABS < TC-TA > >P 1/2  THEN  OOTO  5340 

5330  GOSUB  05300  REM  CALC-LENOTH 

5340  09-AL  TA-TH  :  PA-PH  '  TC-TB  :  PC-PB 

5350  REM  D1 -DISTANCE  FROM  ORIGIN  TO  POINT 

3360  REM  D2-DI8TANCE  FROM  DESTINATION  TO  POINT 

5370  REM  D3-D I  STANCE  FROM  TANKER  BASE  TO  POINT 

5330  REM  09-DISTANCE  FROM  ORIGIN  TO  DESTINATION 

5400  RETURN 

5410  REM 

5420  REM 

5430  REM 

5500  REM  CALC-LENOTH  ALCTHI, PHI, TH2.PH2) 

5600  REM  ALPHA-ANGLE  BETWEEN  THE  TWO  POINTS 
5700  REM  K1 -COSINE  OF  ALPHA 
5800  REM  K3-TANGENT  OF  ALPHA 

5900  Sl-SINCPO  S2-SIN(PR>  D-TC-TA  CD-COSCD)  Cl-COSCPC)  C2-COSCF 
6000  K1=S1*S2*CD+C1»C2  K2-KIT2  :  K5-SQRC1.-K2) 

6100  K3-K5/K1  ■  AR-ATNCK3)  AL-AR4R 

6150  REM  AL-ACTUAL  LENOTHC DISTANCE  BETWEEN  THE  TWO  POINTS) 

6200  RETURN 
6210  REM 
6220  REM 
6230  REM 
6300  REM  RANGE 

6310  REM  CALCULATE  MAXIMUM  RANGE  FOR  EACH  AIRCRAFT 
6320  REM  WHEN  LOADED  WITH  MAXIMUM  FUEL 
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£500  GOSUB  07200  :  REM  RKGMAX.WO) 

€600  Rl-RT  :  GR-OC 

6700  GOSUB  07200  :  REM  RT<GCAP,WO> 

€800  R2-RT  GR-HM 

€900  GOSUB  07200  REM  RT(HMAX,WO> 

7000  R3-0.5*RT 

7010  REM  R3-TANKER  RANGE  TAKING  RETURN  TRIP  INTO  ACCOUNT 

7100  RETURN 
7110  REM 
7120  REM 
7130  REM 

7200  REM  RT<GO,WO> 

7210  REM  CALCULATES  RANGE  FOR  AIRCRAFT  GIVEN  ALLOWABLE  FUEL  AND  CARGO  WEIGHT 

7220  REM  ASSUMES  TRANSPORT  AND  TANKER  HAVE  SAME  PERFORMANCE  CURVES 

7300  RT=' AP+A1*WR+A172*GR>*0R 

7490  RETURN 

7410  REM 

7420  REM 

7430  REM 

7500  REM  SOLVE  SUBPROBLEM  ONE 

7510  REM  FIND  THE  LEAST  COMBINED  FUEL  FOR  TRANSPORT  AND  TANKER 
7520  REM  UNDER  THE  GIVEN  CONDITIONS 
7530  OOSUB  25400  '  REM  FEASIBLE 

7540  REM  UPDATE  DISTANCES  AFTER  SEARCH  IN  SUBPROBLEM  TWO 
76O0  SA=SQR<At2*2#Al*D2>  SB=3GR<AP12+2*A1*D3> 

7700  C=3A+SB  :  UL-WO  :  DL-D1 

7710  REM  C=CONSTANT  USED  IN  LATER  CALCULATIONS 
7300  OOSUB  08200  REM  FNTR<WL,DL> 

7810  REM  FNTR-FUEL  NEEDED  BV  TRANSPORT  TO  TRAVEL  TO  INITIAL  POINT  GIVEN 

7900  LO=FT  HI-GM 

8000  OOSUB  03600  :  REM  SP1 SEARCH (C) 

8010  REM  LOOK  FOR  LEAST  COMBINED  FUEL  COST 
8100  RETURN 
8110  REM 

8200  REM  FNTR<WL.DLXUO,D> 

8210  REM  FNTR-FUEL  NEEDED  BV  TRANSPORT  TO  REACH  REFUEL  POINT 

8220  PEM 
8230  REM 

8300  SC=<AP+Al*WL>f2+2*Al*DL 

8400  SD-SGRCSO  FT=-WL-flP/Al+SDXflI 

8500  RETURN 

8510  REM 

8520  REM 

8530  REM 

8600  REM  SP1 SEARCH <C> 

8610  REM  THIS  IS  A  OOLDEN  SECTION  SEARCH  ALONG  THE  FUEL  RANGE  FROM  JUST  BARE’ 

8620  REM  REACHING  THE  POINT  TO  THE  MAXIMUM  THE  TRANSPORT  CAN  CARRV  AT  TARE OR" 

8700  SE-SQRC5. )  '  TU«<-I . +SE>*.5  AC-0.00010 

8800  G-LO  :  DD-HI 

3900  XA=LO+a.-TU)*<HI-LO> 

9000  XB-LO+TU*<HI-LO> 

9100  TX-XA 

9200  OOSUB  11100  :  REM  TFCKXA.C) 

9300  FA-T1  :  TX-XB 

9400  OOSUB  11100  :  REM  TFCKXB.C) 

9500  FB-T1 

9600  IF  FA>FB  THEN  GOTO  09700 
9650  IF  FA<*FB  THEN  GOTO  10000 

9700  G=XA  XA-XB  XB»G+TU*< DD-O)  ■  FA-FB  TX-XB 
9800  OOSUB  11100  '  REM  TFCICXB.O 
9900  FB-Tl  '  OOTO  10300 

10000  DD-XB  XB-XA  XA-G+C 1 . -TU)B(DD-G)  FB-FA  TX-XA 
10100  OOSUB  11100  REM  TFCKXA.C) 

10200  FA-T1 

10300  IF  ' DD-O-ACX-0  THEN  GOTO  10400  10350  IF(DD-G-AC»0  THEN  GOTO  09600 


10500  SG-SQRCSF)  SH-(C-SG)«+2*A1*D3  SI-SQR<SH> 

10600  H0-SI/A1 -APXA1  F»F8 

10700  IF  FBCFA  THEN  GOTO  11000 

10800  G0-XA  SF«(A+A1*XA>12-2*A1*D1  '  SG»SQR(SF> 

10900  SH=(C-SG>f2+2*Al#D3  SI«SQR<SH>  HO-SI/A1-AP/A1  F-FA 

10910  REM  F-L0UE3T  COMBINED  FUEL  COST  FOR  THE  TRIRL  POINT 

11000  RETURN 

11010  REM 

11020  REM 

11030  REM 

11100  REM  TFCKGO.O 

11110  REM  TFC1 "TOTAL  FUEL  COST  FOR  SUBPROBLEM  ONE 
11120  REM  GIVEN  INITIAL  FUEL  G0  AND  CONSTANT  C 
112O0  SF=<A+Al*TX)t2-2*fll*Dl 

11300  SG-SQR(SF)  :  SH»(C-SO> T2+2*A1*D3  SI-SOR(SH) 

1140O  HO=S I /A 1 -AP/A 1 
11500  Tl-TX+HO 
11510  REM  T1-TFC1 
11600  RETURN 
11610  REM 
11620  REM 
11630  REM 

11700  REM  SOLVE  SUBPROBLEM  TWO 

11710  REM  SEARCHES  FOR  THE  BEST  REFUELLING  LOCATION  GIVEN  INITIAL  FUEL  00 
11720  REM  STARTS  FROM  ORIGINAL  POINT  AND  PROCEEDS  TO  CHECK  DIAGONALS  FOR 
11730  REM  IMPROVEMENT 

11740  REM  CHECKS  EACH  TRIAL  POINT  FOR  FEASIBILITY  BEFORE  EVALUATING 
11750  REM  IN  CASE  OF  TIES  THE  AMOUNT  OF  CHANGE  BETWEEN  POINTS  IS  REDUCED 
11760  REM  UNTIL  THE  TIE  IS  BROKEN  OR  THE  ACCURACY  LIMIT  IS  REACHED. 

1180C  DV=0. 1  :  T9-TH  P9-PH  DZ-0. 1 

11850  TA-TH  PA-PH 

11900  GOSUB  22600  REM  FVSP2 

12000  FF=F2  EP-0.1 

12100  TH-T9+DZ  PH-P9+D2  TV-0. 

12200  GOSUB  25400  REM  FEASIBLE-NOPRNT 

12300  IF  INO50  THEN  GOTO  13400 

12400  GOSUB  22600  ■  REM  FVSP2 

12500  IF  F2CFF  THEN  T9-TH 

12600  IF  F2CFF  THEN  P9«PH 

12700  IF  F2-FF  THEN  TY*1. 

12800  IF  F2CFF  THEN  FF-F2 

12900  IF  fFF=F2>  AND  <TV»0.>  THEN  GOTO  12100 
13000  IF  <FF=F2>  AND  < TV- 1.) THEN  DZ-DZT2 
13100  IF  DZCE2  THEN  GOTO  13300 
13200  IF  <FF=F2>  AND  <TV-1.)  THEN  GOTO  12100 
13300  DZ-EP 

13400  TH-T9-DZ  ;  PH-P9+DZ  ■  TV-0 

13500  GOSUB  25400  REM  FEASIBLE-NOPRNT 

13600  IF  IN 050  THEN  GOTO  14700 

13700  GOSUB  22600  ■  REM  FVSP2 

13800  IF  F2<FF  THEN  T9-TH 

13900  IF  F2<FF  THEN  P9-PH 

14000  IF  F2-FF  THEN  TV-1. 

14100  IF  F2<FF  THEN  FF-F2 

14200  IF  (FF-F2)  AND  (TV-0. >  THEN  GOTO  13400 
14300  IF  (FF-F2)  AND  <TV-1.)  THEN  DZ-DZT2 
14400  IF  DZCE2  THEN  OOTO  14600 
14500  IF  (FF-F2>  AND  <TV-1.>  THEN  OOTO  13400 
14600  DZ-EP 

14700  TH-T9-DZ  PH-P9-D Z  TV-0. 

14800  GOSUB  25400  REM  FEASIBLE-NOPRNT 

14900  IF  INO50  THEN  OOTO  16000 
15000  OOSUB  22600  REM  FVSP2 

15100  IF  F2CFF  THEN  T9-TH 
152O0  IF  F2<FF  THEN  P9-PH 
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15400  IF  F2CFF  THEN  FF-F2 

15500  IF  (FF-F2)  AND  (TV-0.)  THEN  GOTO  14700 
15600  IF  '.FF-F2)  AND  (TV-1.)  THEN  DZ»DZt2 
1570O  IF  D2CE2  THEN  GOTO  15900 
15300  IF  (FF-F2)  AND  <TV-1.>  THEN  OOTO  14700 
15900  DZ=EP 

16000  TH-T9+DZ  PH-P9-DZ  TY-0. 

1610O  GOSUB  25400  ■  REM  FEAS I BLE-NOPPNT 

162O0  IF  INO50  THEN  GOTO  17300 

16300  GOSUB  22600  :  REM  FVSP2 

1640O  IF  F26FF  THEN  T9-TH 

165O0  IF  F2<FF  THEN  P9-PH 

16600  IF  F2-FF  THEN  TV-1. 

16700  IF  F2CFF  THEN  FF-F2 

16800  IF  (FF-F2)  AND  (TV-0.)  THEN  GOTO  16000 
16900  IF  cFF»F2)  AND  (TV-1. >  THEN  DZ-DZt2 
17000  IF  DZCE2  THEN  OOTO  17200 
17100  IF  (FF-F2)  AND  (TV-1.1  THEN  GOTO  16000 
17200  D2-EP 

17300  IF  EPCE2  THEN  OOTO  17550 
17400  EP-EP4DV  D2-EP 
17530  GOTO  12100 
1 7550  TH=T9  PH-P9 

17560  GOSUB  50500  REM  CONVERT  SOLUTION  TO  LONO  AND  LAT  TO  SHOW  ITERATION? 
17570  PRINT  “SUBPROBLEM  TUO  TRIAL  SOLUTION" 

17575  PRINT  L9; D7*. H8, D8»- F201000 

17590  REM  PRINT  SOMETHING  TO  SHOW  OPERATOR  THAT  THE  PROGRAM  IS  PROGRESSING 

17600  RETURN 

17610  PEM 

17620  PEM 

17630  REM 

1S00O  REM  SOLVE  SPECIAL  CASE 

18010  PEM  REFUELLING  IS  NOT  NEEDED  HERE. 

13020  REM  THEREFORE  THE  FUEL  NEEDED  FUNCTION  IS  SOLVED  AND  LOCATION  FOP 

13030  PEM  REFUELLING  IS  SET  TO  0 

18100  SC«(AP+A1*W0>12*2*A1*09 

18200  SD-SQR(SC)  G0=-WO-AP/A1*:3D/M1 

18300  L9=L2  D7*=3L*  H8-H2  D8B-SM*  HO-0  FT-8  F-00  K-l 

13400  RETURN 

13410  PEM 

18420  PEM 

18430  PEM 

19000  PEM  SPECIAL  CASE  WHERE  DISTANCE  I  PI 72 
19100  AM-TC  PM-PC  AN-TA  :  PN-PA 
19200  IF  TC-TAC0  THEN  TC-AM+( ABS( AM-AN>)/2 
1930O  IF  TC-TA30  THEN  TC-AM-CAM-AN V2 
19400  PC=(PM+PNV2 

19500  GOSUB  20000  PEM  CALC-LENGTH 

19600  LQ=AG>  TA-TC  PA-PC  TC-AM  PC-PM 

197O0  GOSUB  20000  PEM  CALC-LENGTH 

19800  LP=AQ  AL-LQN.P  TC-AM  PC=PM  TA-AN  PA-PN 

19900  RETURN 

19910  REM 

19920  PEM 

19930  REM 

20000  REM  CALC-LENGTH  FOR  DISTANCE  >  PI/2  ROUTINE 

20100  S1=SIN(PC>  ■  S2-SINCPA)  D-TC-TA  :  CD-COS<D)  Cl-COSCPC)  C2-C0S'?A 
20200  K 1 =S 1 *S2*CD+C 1 #C2  K2-K1T2  K5-SQRC 1 . -K2) 

20300  K3-KS/K1  AR-ATN(K3>  AO-ARPR 

20400  RETURN 

20410  REM 

20420  REM 

20430  REM 

22600  PEM  F VSP2  <  THSP2 . PHSP2  > 

22610  REM  FVSP2-FUEL  COST  FUNCTION  FOR  SUBPROBLEM  TWO 

''-•6.70  RFM  I.IHTPM  Tg  PI  TnuTI  V  HTFFCPCNT  FROM  81  IRgBOBLEH. JJk g-  JECaUSE. XlE - 


22630  REM  THE  CONSTRAINTS  BEING  RELAXED. 

22700  GOSUB  04400  REM  DISTANCE 

22710  REM  NEW  DISTANCES  MUST  BE  CALCULATED  EVERV  TIME  THE  REFUELLING  POIHT  ' 
22720  REM  CHANGED. 

72900  GOSUB  24600  '■  REM  FCTR<THSP2> 

27310  REM  FCTR-FUEL  COST  FOR  THE  TRANSPORT-AGAING  DIFFERENT  BECAUSE  OF  PEL-‘r 
32820  REM  CONSTRAINTS. 

22900  WP-140+ ( G0-FX )  WL-0.0  DL-D3 

23000  GOSIJB  25000  REM  FNTA<0,D3> 

23050  REM  FUEL  NEEDED  IV  TANKER  TO  RETURN  AFTER  REFUELLING 

23100  F3=FV 

23200  ML  =140  DL»D 2 

23300  GOSUB  08200  REM  FNTRO40.D2) 

2335©  REM  FUEL  NEEDED  IV  TRANSPORT  TO  GET  TO  DESTINATION  AFTER  REFUELLING 

23400  F4=FT 

23500  W=F3+F4-<G0-FX) 

23550  REM  NET  REQUIREMENTS  OF  TANKER  TO  CARRV  AS  CAROO 

23600  UL=MP  DL=D1 

23700  GOSUB  08200  REM  FNTR<I4P,D1> 

23750  REM  FUEL  NEEDED  BV  TRANSPORT  TO  GET  TO  POINT  WITH  CARGO  AND  UNUSED  F'JE 

23800  F5=FT  14L-UI  :  DL=D3 

24100  OOSUB  25000  REM  FNTA(W,D3> 

24110  REM  FUEL  NEEDED  BV  TANKER  TO  GET  TO  REFUELLING  POINT  WITH  NET  FUEL -I  H 

24200  F7=FV 

34400  F2-F5+F4+F7+F3 

24410  REM  F2-FUEL  COST  FOR  SUBPROBLEM  TWO 

24500  RETURN 

24510  REM 

24520  REM 

24530  REM 

24600  REM  FCTR(THSP2) 

24610  REM  FCTR-FUEL  COST  FOR  TRANSPORT 

24700  SV=<AP+A1*<WO*G0>>  :  SW»SVT2-2*A1*D1 

24800  SX=SQR<SW>  EX-O0+WO+ AP7A1 -SX7A 1 

24900  RETURN 

24910  REM 

24920  REM 

24930  REM 

25000  REM  FNTA<WL.DLXWO,D> 

25010  REM  FNTA=FUEL  NEEDED  BV  THE  TANKER 

25020  REM  ASSUMES  TANKER  AND  TRANSPORT  HAVE  IDENTICAL  PERFORMANCE  CURVES. 

25100  SC=<AF'+H1*WL>T2+2»A1*DL 

25200  SD-SQR<SC>  FV—UL-AP/A1+SD/A1 

25300  RETURN 

25310  REM 

25320  REM 

25330  REM 

25400  REM  FEA3 1 BLE-N0PR1NT 

25410  REM  SAME  AS  PREVIOUS  ROUTINE  BUT  DOES  NOT  PRINT  ANV  MESSAGES. 

25500  TA=TH  PA=PH 

25600  GOSUB  04400  REM  DISTANCE 

25700  GOSUB  06300  REM  RANGE 

25710  WL=WO  DL=D1  OOSUB  08200  REM  FNTRCWO . DL ) 

25800  IN=0 

25900  IF  ■  Dl>Rl  >  OP  <FT>GM’>  THEN  IN-10 
26000  IF  D2>P2  THEN  IN=20 
26100  IF  D3>R3  THEN  IN=30 

26200  IF  I  HO  10  AND  INO20  AND  INO30  THEN  IN-50 

26300  RETURN 

26310  REM 

26320  REM 

26330  REM 

2S00O  REM  CONVERSION  TO  PAD IANS  FOR  COORDINATES 

23010  REM  ALL  COORDINATES  ARE  DEFINED  IN  RELATION  TO  THE  ORIGIN. 

23020  REM  THIS  ALLOWS  THE  CATEGORIZATION  OF  MOST  PROBLEMS. 
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23040  REM  INTERNATIONAL  DATE  LINE  OR  ACROSS  THE  EQUATOR.  ALSO  IF  ROUTES  ALO! 
28050  REM  THE  SOUTHERN  LATITUDES  ARE  TO  BE  CONSIDERED. 

23100  L7=L1-L1 

28200  IF  <SJ*=''W“)  AND  <SL*="U")  AND  <SN#="U")  AND  <SQ*»*U">  THEN  L4«L1-L2 

23300  IF  <SJ*=“U”>  AND  <SL*»"W")  AND  <SN*="W">  AND  (SQta'U")  THEN  L5-L1-L3 

28400  IF  <SJ*=“W>  AND  <SL*»"W">  AND  <SN*«"W>  AND  <SQ*»"U">  THEN  L6-U-CA 

28500  IF  <SJ*«”M“>  AND  <SL*="UH)  AND  <SN*«"W)  AND  (SOM-U")  THEN  02*1 

28600  IF  <02=1 >  AND  <L1>L2)  AND  CL1>L3)  AND  <LI>CA>  THEN  QQ-1 
23650  IF  <0Q=1)  THEN  GOTO  50000 

28200  IF  <SJ*="U")  AND  <SL*="E")  AND  <SN*«"W“)  AND  (SQ«-"U">  THEN  L4-L1+L2 

28800  IF  <SJ*="U")  AND  <SLf=“E")  AND  <SNS»''W")  AND  <SQ*-“U“>  THEN  L5-L1-L9 

28800  IF  <SJ*»"W*>  AND  <SL*«"E")  AND  <SN**"W')  AND  <SQ*a*U“>  THEN  L6=L1-CA 

28000  IF  <SJ*="W“>  AND  <SL*=“E")  AND  <SN*=“U">  AND  <SQS»"M")  THEN  Q2=2 

29050  IF  <02=2)  AND  <L4<=180>  AND  <L1>L3>  AND  CLOCA)  THEN  OQ-2 
29075  IF  <Q0=2>  THEN  GOTO  50000 

30000  IF  )  AND  <SL*=“E“>  AND  <SN*=“U")  AND  <SQ*«“E“>  THEN  L4=L1+L2 

3010O  IF  <3J*»“W“>  AND  <3LM*E">  AND  <SN*="U")  AND  (SQ*«"E">  THEN  L5=L1-L3 

3O200  IF  <SJ*=“W)  AND  <SL*«*’E")  AND  <SN$«"U”>  AND  <SQ*»"E")  THEN  L6=L1+CA 

30300  IF  <S.T*=“M")  AND  <SL**“E‘'>  AND  <SN*="W")  AND  <8Q*="E")  THEN  QZ=3 

30400  IF  < 02=3  >  AND  <L4<=I80>  AND  <L1>L3)  AND  <L6<«19e>  THEN  QQ-3 
30450  IF  <Q0=3>  THEN  GOTO  50000 

30500  IF  <SJ*="W")  AND  <SL*="E">  AND  <SN*="E">  AND  <SQ#»"14")  THEN  L4=L1+L2 

30600  IF  <SJ*="U")  AND  <SL*="E">  AND  <SN*»"E">  AND  <SQ*="U">  THEN  L5-L1+L3 

30700  IF  <3J*="W')  AND  <SL*="E")  AND  <SN*»"E")  AND  <8Q*»"W">  THEN  L6«L1-CA 

3080O  IF  <SJ*="W")  AND  <SL*="E“)  AND  <SN*="E")  AND  <SQ*«"W“>  THEN  QZ=4 

30900  IF  <QZ«4>  AND  <L4<=180)  AND  <L5<=180)  AND  <CACL1>  THEN  QQ-4 
30950  IF  <00=4  >  THEN  GOTO  50000 

31000  IF  <SJ*="U")  AND  <SL*= "E” )  AND  <SN$«”E“)  AND  CSQ*=“£">  THEN  L4=L1+L2 

31100  IF  .-SJ*="W")  AND  <SL*«“E">  AND  <SNt»"E“)  AND  CSQ*»"E">  THEN  L5«L1+L3 

31200  IF  <SJf ="W" )  AND  <SL*="E”)  AND  <SN*»"E">  AND  <SQ#-"E">  THEN  L6-L1+CA 

31300  IF  ■’SJ*="W")  AND  <SL*="E")  AND  <SN$»»E")  AND  <SQ*«"E">  THEN  Q2=5 

31400  IF  <0Z=5>  AND  <L4<«I80)  AND  <L5<»180>  AND  <L6<-180>  THEN  00=5 
31450  IF  <00=5 >  THEN  GOTO  50000 

31500  IF  <SJ*=''E“)  AND  <SL*="E")  AND  <SN*="E">  AND  <SQ*»“E">  THEN  L4-L2-L1 

31600  IF  <SJ*="E")  AND  <SL*="E")  AND  <SN*»"E'’>  AND  <SQ*«"E'’>  THEN  L5=L3-L1 

31700  IF  <SJt=”E" )  AND  <SLf»"E")  AND  <SN*="E")  AND  <SQ*«"En>  THEN  L6=CA-L1 

31300  IF  <SJ*="E")  AND  <SL*="E”)  AND  <SN*=”E*)  AND  CSQS-'E”)  THEN  02=6 

31900  IF  <0Z»6>  AND  <L2>L1)  AND  <L3>L1)  AND  <CR>L1>  THEN  00-6 
31950  IF  <0Q=6>  THEN  GOTO  50000 

33000  IF  <SJ*="W")  AND  <SL*«"U")  AND  <SN*="W">  AND  <SQ*»"U")  THEN  L4=L2-L1 

3310O  IF  <SJ*="W")  AND  CSL*="U")  AND  <SN*="U">  AND  <SQ*»"U")  THEN  L5-L3-L1 

33200  IF  <SJ»="W">  AND  <SL*="U")  AND  <SN*="U" )  AND  <SQ*=“W">  THEN  L6=CA-L 1 

33300  IF  <SJ*="W")  AND  <SL*="U">  AND  <SN*="I4")  AND  <SQ*=“W")  THEN  QZ=7 

33400  IF  <02=7)  AND  <L2>L1)  AND  <L3>L1>  AND  <CA>L1>  THEN  QQ«7 
33450  IF  <00=7 >  THEN  GOTO  50000 

33500  IF  <SJ$=“E")  AND  <SL*=”W")  AND  <SN»«"U">  AND  <SQ*=“W")  THEN  L4=L1+L2 

33600  IF  <SJ*="E")  AND  <SL*="U">  AND  <SN*="W“)  AND  <SQ#="U")  THEN  L3=L1+L3 

33700  IF  <SJ*="E" >  AND  <SLf=''U">  AND  <SN*="U">  AND  <SQ*="U'')  THEN  L6=L1*CA 

33800  IF  <SJ#»"E")  AND  <SL*="W">  AND  <SN#="H">  AND  CSQ*=*W">  THEN  Q2«S 

33900  IF  <02=8)  AND  <L4O180)  AND  <L5<=130)  AND  <L6O=180)  THEN  00=8 
33950  IF  <00=8)  THEN  GOTO  50000 

34000  IF  <3Jt*"E" )  AND  <SL*="W”)  AND  <SN*«"W")  AND  <SQ*="E">  THEN  L4-L1+L. 

34100  IF  <SJ*="E")  AND  <SL*="W")  AND  <SN*-”M">  AND  <SQ*«"E")  THEN  L3«L1*L: 

34200  IF  <SJt="E")  AND  <SLM“W">  AND  <SNt="H")  AND  <8Qf="En>  THEN  L6«L1-CA 

34300  IF  <SJ*=“E")  AND  <SL*=’’W")  AND  <SN*="U”)  AND  <SQ*=“E“)  THEN  Q2=9 

34400  IF  <QZ=9)  AND  <L4<=180)  AND  <L5<=180)  AND  <CA<L1)  THEN  00=9 
34450  IF  <QQ=9)  THEN  GOTO  50000 

34500  IF  <SJf="E")  AND  <SL*«"tJ")  AND  <SN*="E">  AND  <SQ*=”U")  THEN  L4=L1*L2 

34600  IF  <SJ*="E” >  AND  <SL*«"U" )  AND  <SN*="E")  AND  <SQ*»"H")  THEN  L5=L1-L3 

34700  IF  'SJ»="E">  AND  <SL*="U")  AND  <SN*="E")  AND  <SQ»="H">  THEN  L6=L1*CA 

24800  IF  <SJ*="E“)  AND  <SL»=“W")  AND  <SN»="E")  AND  <8Q*="U”)  THEN  Q2=10 

34900  IF  <02=10)  AND  <L4<=180>  AND  <L3<l.l)  AND  <L6<«180)  THEN  00=10 
34950  IF  <00=10)  THEN  GOTO  50000 

35000  IF  ''SJ*»"E")  AND  <SL*="W">  AND  <SN*»"E")  AND  <SQt»"E")  THEN  L4=L1+L- 

35100  IF  <SJ*="E")  AND  <SL*="W")  AND  <SN*»"E">  AND  <SQ#=”E’'>  THEN  L5=L1-L: 
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35300  IF 
35400  IF 
35450  IF 
35500  IF 
35600  IF 
3570O  IF 
35800  IF 
35900  IF 
35950  IF 


'.SJ*«"E“>  AND  <SL**"U")  AND  <SN»-"E“)  AND  <8Q*«"E">  THEN 
<Q2=1 1 >  AND  CL4<=180)  AMP  <L3<L1)  AND  <CA<Lt>  THEN  QQ-11 
QQ* 1 1 >  THEN  OOTO  50000 


<SJ*«"E")  AND  <SLB*''E")  AND  <SN*-"E")  AND  <SQ«»"E">  THEN 

:SJ*»"E")  AND  CSL*"“E")  AND  <SN**"E“>  AMD  <SO0*“E">  THEN 

<SJ**"E">  AND  <SL*-"E")  AND  <SN*»"EM>  AND  (SQ*»“E">  THEN 

|,SJ*«"E" )  AND  <SLt*"E")  AND  <SN*-”E“)  AND  <3Q*-"E">  TWN 

<Q2*  12)  AND  <L2<L1>  AND  (L3<Lt)  AND  <CACL1>  THEN  00*12 
..Q0*12)  THEN  OOTO  50000 


50000  T0=L?  :  PO*P I /2-H 1 *P 1/ 1 80 


50100  TD=L4*PI/180  :  PD-PI/2-H2BPI/180 
50200  TE=L5*PI/180  :  PS*PI/2-H3*PI/180 
50300  TH=L6*PI/180  PH*PI/2-CB*PI/180 

50310  REM  CONVERT  THE  COORDINATES  TO  RADIANS 

50320  PEM 
50330  PEM 
50340  REM 


02*11 


L4*L1-L2 

LS-L1-L3 

L6=L1-CP 

Q2-12 


50350  RETURN 

50500  PEM  CONVERT  SOLUTION  TO  LONOITUDE  AND  LATITUDE 

50510  REM  THIS  ROUTINE  WOULD  ALSO  NEED  TO  BE  EXPANDED  IF  MORE  CASES  ARE  TO  I" 
50520  PEM  CONSIDERED. 

50600  H8=90-PH* 180/FI  DeB-'N" 

5O700  IF  <Q0<*5 1  OR  (<0Q>-8)  AND  •0Q<«12))  THEN  L9»L1-<TH*1 80/PI > 

50800  IF  090=6)  OR  <Q0*7)  THEN  L9*L1 ♦<TH*I 80/PI ) 

50900  IF  <Q0O5>  AND  <L9>*0)  THEN  D7B»"U" 

51000  IF  <0GO5)  AND  <L9C0)  THEN  D?B-"E" 

51100  IF  '  QQ<*5)  AND  <L9<0)  THEN  L9-ABSIL9 ' 

51200  IF  00=6  THEN  D7B-“E" 

51300  IF  00=7  THEN  D7**"W" 

51400  IF  < (00>*8)  AND  CQ0<*11))  AND  <L9>*0>  THEN  D?B*"E" 

515O0  IF  «0Q>-8>  AND  f0QC*It>)  AND  CL9C0)  THEN  D7B-’H” 

51600  IF  ^QQ>=8)  AND  'WX-ll))  rnd  <l9<0>  THEN  L9-ABSCL9) 

51700  IF  Q0*12  THEN  D7B*“E" 

51800  RETURN 
51810  REM 
51320  PEM 
51830  REM 

60500  REM  PRINT  MAP  OF  U.S.  AND  EUROPE 

60501  PRINT  "d" 

60507  P0KE1188.43  POKEl 191 . 43 : POKEl 194. 43  POKE 1 197 , 43  POKE  1290. 43  POKE 1284 

60508  POKE  1285 . 43  POKEl 286 • 43 ■ POKEl 21 2. 43 1 POKEl 2 13. 43  POKEl 227, 43  P0KE122S  4: 

60509  POKE  1231 , 43  POKE  1232, 43  POKE  1233. 43  P0KE1234.43  P0KE1238.43  POKE1240  4; 

60510  POKE  1251,43  POKE  1 254 ,43  POKE  1 255 , 43  P0KE1267.43  P0KE1271.43  P0KE127^  J 

60511  POKE 1278. 43  POKE 1279, 43  P0KE129I ,43  P0KE1294.43  P0KE1293.43  POKE130-  l\ 

60512  POKE 131 1,43  POKE  131 3, 43  P0KE1314,43  P0KE1315, *3  POKE 1320, 43  POKEl 329  4T 

60513  POKE  1 332, 43 : POKE  1 333 , 43 : POKE  1 334 , 43  POKE 1333, 43  POKE 1340, 43  POKEl 34 Q  4 

60514  POKE  1331 .43  POKE  1 355 , 43  POKE  1368. *3  POKE  1370, 43  POKEl 372. 43  POKE  1 37?  ■: 
60513  POKE 1374, 43  POKE  1390. 43  P0KE1391 - 43  P0KE1396, 43  POKE1407.43  POKE1410  4: 

60516  POKE  141 1 ,43  POKE  1430- 43  P0KE1436, 43  P0KE1448, 43  P0KE1449, 43  POKE1450  >' 

60517  POKEl 477 , 43  POKE 1489 ■ 43 : POKE  15 13, 43  P0KE151343  P0KE1314.43  POKE1516  i: 

605 1 5  POKE  1 528 , 43  POKE  1 333 . 43  POKE 1567,43  POKEl 373 , *3 : POKE 1374 , 43  POKE 1374  4 ' 

60519  POKE  1575, 43 : POKE  1592, 43 : POKE  1603, 43  POKE 1607,43  POKE 1608, 43  POKE 1609  4? 

60520  POKE  1610, 43 : POKEl 61 1,43’ POKE  161 2, 43  P0KE1614, 43  P0KE1613, 43  P0KE1616  4 

6052 1  POKE 1617,43' POKE 1618,43: POKEl 631 , 43 : P0KE1648 , 43 : POKEl 653, 43 : POKE I 635  4 ; 

60522  POKE 1 658 , 43  POKE 1 670 , 43 : POKE 1687, 43 : POKE 1 694 , 43 ; POKE 1696.43: POKE 1 69”  4 : 

60523  POKE  1 698 , 43 : POKE 1 702 , 43  POKE 1 706 , 43  POKE 1707,43: POKE 1708,43: POKE 170?  4 ; 

60524  POKE  1710,43  POKE 1727,43  ’ POKE 1738,43: POKE 1 739, 43 : POKE 1 742, 43 : POKE  1742. 4' 

60525  POKE  1 743 , 43  POKE  1 744 , 4  3  POKE  1 746 , 43 :  POKE  1730,43:  POKE 1731, 43  POKE  1 766  4  ;■ 
60326  POKE  1779, 43 : POKE  1785  43  POKE 1787, 43 : POKEl 788, 43  POKE 1789, 43  POKE 1791 .4? 

6052 7  POKE 1 803 • 43 : POKE 1 8 1 9 , 43 : POKE 1 820 , 43 : POKE  1826,43: POKE  1827,43: POKE  1 828  4 

60528  POKE 1 829 , 43  POKE 1 830 , 43  POKE 1843 , 43 : POKE 1060.43: POKE  1861,43: POKE  1 862  ’ ’ 

60529  POKE  1863, 43 : POKE  1 869 . 43 ■ POKE 1 870, 43 : P0KE1872, 43 : POKE  1073, 43 : POKE  1874.4 ‘ 

60530  POKE  1886, 43  POKE  190 1,43  POKE1902,43  POKE1910,43:POKE19U,43  POKE1916  J 

60531  POKE  1926, 43 1  POKE  1927, 43 : POKEl 928, 43 : POKEl 929, 43  POKE 1938, 43 : POKE 1942  4 : 

60532  POKE  1931,43: POKE 1937, 43  P0KE1938. 43  POKE 1971, 43 : POKE 1972, 43 : POKE 1981  - 
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-..O'. ?4  FOR-  U=1992T01999  POKEU, 102  NEXTUFOR  0201 3TO2820 : POKEU, 162 : NEXTU 
80525  FOP  U=  1 932T0 1 956  POKEU .  1 02 : NEXTU ; FOR  0 1973T01 960 : POKEU.  102  NEXTU 
60536  FOP  U-1912T01913  POKEU,  102 : NEXTU : FOR  U-1931T01941 : POKEU, 182 :NEXTU 
-'0337  FOR  U-1886TO1900  POKEU,  102  NEXTU  FOR  0 1 846T0 1 039 ■  POKEU,  182 NEXTU 

60338  fop  IJ-1006TO1810  POKEU,  132  NEXTU  FOR  01821  TO 1623  POKEU,  102 : NEXTU 

8033?  FOR  U" 1767T01 778 ■ POKEU, 102 : NEXTU : FOR  O 1 780TO 1 783  POKEU .  182 : HEXTU 
8.0340  FOR  IO1728T01737  POKEU,  102  HEXTU  FOR  0 1 664T0 1669 : POKEU,  182 : 1CXTU 
80341  FOP  U-1S88T01S93  POKEU, 102  NEXTU -FOR  O 1699T01 70 1  POKEU, 182 : NEXTU 
8034?  FOP  O1624TO1630  POKEU, 102  NEXTU- FOR  O 1 649T0 1 632 • POKEU ,182 HEXTU 
8.0543  FOP  O1659T01663  POKEU,  102  NEXTU  FOR  O1504TO1391  POKEU,  182:NEXTU 
80344  FOR  U-1619T01623  POKEU, 102  NEXTUFOR  O1344T01332  POKEU. 182:NEXTU 
8.0343  FOR  U=1376T01383 '  POKEU,  102  NEXTU  FOR  U«1304TO1312: POKEU,  102 :  NEXTU 
8.0346  FOP  U=1529T01543 ; POKEU,  102 : NEXTU  FOR  U-1464T01476 : POKEU,  182 'NEXTU 

80547  FOR  U=1490TO 1503  POKEU, 102  NEXTU  FOR  U-1424T01429  POKEU. 182  NEXTU 

80548  FOP  U«1431T01433  POKEU, 102  NEXTU  FOR  U-1431T01463 : POKEU. 182  NEXTU 

80349  FOP  IJ-1384T01389  POKEU  ,102-  NEXTU :  FOR  U-  1392T01 393  POKEU,  182  NEXTU 

60330  FOP  U*1412T01423  POKEU,  102  NEXTU  FOR  U-1373T0 1383  POKEU. IK^NEXTU 

80551  FOP  U»1336T01 343  POKEU,  1 02 : NEXTU : FOR  U=1296TO1303POKEU.182NEXTU 

80552  FOR  U-1256T01263  POKEU, 102  NEXTU  FOR  U-1216T01223  POKEU. 182: NEXTU 

80553  FOR  U«1184T01187  POKEU, 102  NEXTU  FOR  U» 1 224T0 1 226  POKEU , 1 02: NEXTU 

60554  FOP  U=1264T01266  POKEU, 102 : NEXTU : FOR  U=1384T01386:  POKEU. 182 : NEXTU 
60353  FOR  U=1344T01347  POKEU, 102  NEXTUFOR  U-1332T0 1334  POKEU, 182NEXTU 
80536  FOR  O1192T0I193  POKEU,  102  NEXTU  FOR  U-1198T01 199  POKEU, 182  NEXTU 
60337  FOP  U*1213T01214  POKEU, 1 02 : NEXTU : FOR  U“1292T01293: POKEU, 102 : NEXTU 
60353  FOP  U=1408TO1409 : POKEU, 102 ; NEXTU : FOR  U-1368T01369  POKEU, 102  NEXTU 
6.0559  FOP  U=1704TO1705  POKEU,  1 02 :  NEXTU :  FOR  U-  1740TO1 74 1  :  POKEU,  102  NEXTU 
60560  POKE 1239 ,102  P0KE1252, 102 : P0KE1233, 102 ■ P0KE1272, 102 ■ P0KE1312, 102 
6036 1  POKE 1 369 , 1 02  POKE 1703, 102' POKE 1 745 , 1 02  POKE! 786 , 1 02 

80600  7=24- 1 NT <Hl/3. 55 

80603  IF  SJ»=“ir  THEN  H«27-INT(Ll/4) 

60610  IP  S.T*="E"  THEN  H*27+INTCL174) 

60615  POKE  1024+H+404V, 15 
60620  POKE  53296+H+40RV , 0 
60625  7=24- 1NT (H2/3. 5) 

60630  IF  $L*="U"  THEN  H»27-INT<L2/4) 

60635  IF  ?L*="E"  THEN  H=27*INT<L2/4) 

60640  POKE  1024+H+40*V. 4 


60643  POKE  55296+H+40*V,O 
60630  7=24- 1 NT  8H3/3. 5  > 

60653  IF  SMT="M"  THEN  H=27-INT(L3/4) 

60660  IF  SN*="E“  THEN  H=27+INT<-L3K4> 

60663  POKE  1 O24+H+404V , 2 
60670  PONE  332?6+H+40*V,0 
60673  V=24-INT<H873.3'> 

60680  IF  B7*="W"  THEN  H«27-lNT<L?/4) 

60685  IE  P7f =“E"  THEN  H«27+INTa?743 

60690  POKE  1024+H+404V, 18 
80695  POKE  55296+H+404V, 1 

60700  PEN  PLACE  COORDINATES  ON  NOP  FOP  OR  IOIN,  DEST  TANKER  BASE .  ANP  SOU 
60710  PEN  POINT  FOP  REFUELLING 
60720  8 TOP 
60730  RETURN 


81000  STOP 


82000  ENP 


PEAPV. 


0  .V 


KEY  FOR  TRIAL  RUNS  CHART 


t 


LOCATION 

LONGITUDE 

LATITUDE 

New  Jersey 

75W 

40N 

Delaware 

75W 

38N 

North  Carolina 

78W 

35N 

Puerto  Rico 

66W 

18N 

Azores  Islands 

25W 

37N 

Iceland 

20W 

65N 

Germany 

10E 

50N 

Turkey 

30  E 

40N 

Saudi  Arabia 

47E 

25N 

Egypt 

28E 

30  N 

England 

OE/W 

52N 

WEIGHT 

CODE 

100,000 

lbs 

1 

200,000 

lbs 

2 

FUEL  ACCURACY 

100  means  to  within  100  lbs 


POSITION  ACCURACY 
1  means  to  within  1  degree 


TRIAL  RUNS 
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